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Neutrino oscillations have been established in solar and atmospheric neutrinos, but a third signal from the LSND 
experiment is incompatible with three Standard Model neutrinos. The MiniBooNE experiment can confirm or 
refute the LSND oscillation signal with 1 x 10 21 protons on target. While working towards the oscillation result, 
MiniBooNE will accumulate more than 1 x 10 6 neutrino interactions in the to 2 GeV range which will greatly 
increase the world's knowledge of neutrino cross sections in this energy regime. Preliminary results on the 
MiniBooNE charged current quasi-elastic analysis are presented and compared to the NUANCE, NEUT, and 
NEUGEN cross section Monte Carlos. 



1. Introduction 

Neutrino masses have been conclusively estab- 
lished via positive neutrino oscillation measure- 
ments of solar and atmospheric neutrinos £Q, |2J. 
A third and unconfirmed positive result comes 
from the LSND experiment, a short baseline ac- 
celerator 77 M — > ]7 e oscillation search [jJJ. The 
MiniBooNE experiment at Fermi National Ac- 
celerator Laboratory will confirm or refute the 
LSND result with higher statistics and different 
sources of systematic error. 

MiniBooNE will accumulate more than 1 x 10 6 
neutrino interactions on a CHi target. There 
is a dearth of neutrino cross section data on 
heavy targets in the MiniBooNE energy range, 
< E v < 2 GeV. Further, this is a com- 
plex region theoretically since both charged cur- 
rent quasi-elastic (CCQE) and resonance pro- 
cesses contribute in roughly equal proportions. 

The description of this region differs signifi- 
cantly among neutrino cross section Monte Car- 
los. Collaboration between the authors of the 
NUANCE g|, NEUT0, and NEUGEN[H| Monte 
Carlos and the MiniBooNE cross sections working 
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group has made detailed comparisons of CCQE 
event rates and kinematics possible for the first 
time. Confronting the predictions of these cross 
section Monte Carlos with the MiniBooNE CCQE 
data set provides valuable input to the develop- 
ment of simulations of this complex regime. Pre- 
liminary indications point to deficiencies in the 
modeling of low Q 2 interactions common to all 
cross section Monte Carlos considered here. 

2. MiniBooNE Overview 

The MiniBooNE neutrino beam is produced 
from 8.89 GeV/c protons incident on a thick 
beryllium target located inside a magnetic fo- 
cusing horn. Typical operating conditions are 
4 x fO 12 protons per pulse, at 3 - 4 Hz with a 
beam uptime of ~ 88%. The beam spill duration 
is f.6 /is. The focusing horn increases the neu- 
trino flux at the detector by a factor of ~5, and 
can operate in both positive and negative polar- 
ities for v and V running; currently positive sign 
mesons are selected. Mesons produced in the tar- 
get decay in a 50 m long decay pipe. The neutrino 
beam is primarily composed of from ir + decay 
with a 0.5% contamination from v e . The mean 
neutrino energy is ~ 700 MeV; the neutrino flux 
is shown in figure ^ 

The MiniBooNE detector is a 6.f m radius 
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Figure 1. Preliminary MiniBooNE neutrino flux 
Monte Carlo prediction vs. E v (GeV). 

sphere filled with mineral oil (CH2). There 
are 1280 inward-facing "tank" PMTs, and 240 
outward-facing "veto" PMTs. Particle identifi- 
cation depends upon both prompt Cerenkov and 
time-delayed scintillation light. Neutrino induced 
events are identified by requiring that the event 
occur within the beam spill, have fewer than 6 
veto PMT hits, and have greater than 200 tank 
PMT hits. With these simple cuts the cosmic ray 
background is reduced to less than 0.1% of the 
beam- induced neutrino signal 2 . A fiducial vol- 
ume cut at R < 5 m is also typically required to 
ensure good energy reconstruction. 

3. Event Rate Prediction 

The event rate prediction is based on the prod- 
uct of neutrino flux and cross section. The Mini- 
BooNE neutrino flux is primarily produced from 
7r + decay in flight. Therefore a detailed under- 
standing of the 7r + production in p — Be collisions 
is necessary. The neutrino cross section at Mini- 
BooNE energies has contributions from a number 
of different processes jS|. Disentangling the var- 
ious contributions to the total cross section is of 
theoretical interest in this energy regime, and is 
important to oscillation measurements. 



2 The cosmic ray rejection is demonstrated in reference |7|. 



3.1. Flux Prediction 

At MiniBooNE, the relevant ranges of ir + pro- 
duction momenta and angles are 1 < p^ < 
4 GeV/c and < 9 n < 0.2 radians respec- 
tively. There are no existing measurements at 
8.89 GeV/c proton momentum. However there 
is a limited amount of relevant production data 
from past experiments [H] at 10, 12, and 19 
GeV/c. To address the paucity of production 
information, MiniBooNE collaborators have ana- 
lyzed data at 6, 12, and 17 GeV/c from the BNL 
E910 experiment |10j . In addition, 20 million 
triggers were collected with a replica MiniBooNE 
target and an 8.89 GeV jc proton beam at the 
CERN HARP experiment 11 . The analysis of 
this data is currently in progress, and the result- 
ing p Be — > 7r + X cross section measurement 
will be used for the final MiniBooNE flux pre- 
diction. Currently, the neutrino flux is modeled 
with a GEANT4 based Monte Carlo [T2] and an 
external parameterization of the p Be — > n + X 
cross section. The parametrization comes from a 
global fit of existing production data in the range 
10 < Pproton < ITGeV/c to the Sanford-Wang 
model |13) . 

3.2. Cross Section Prediction 

The NUANCE Monte Carlo is used to predict 
the neutrino interaction cross sections. At Mini- 
BooNE neutrino energies the cross section has 
contributions from charged current quasi-elastic 
scattering (39% of the total event rate), charged 
current resonance production (25%), neutral cur- 
rent elastic scattering (7%), and neutral current 
7T° production (7%). For the —>> v e oscillation 
analysis, the most important processes are CCQE 
scattering which affords a precise measurement of 
the neutrino energy, NC 7r° production which is a 
large background to a v e signal |7j, and NC elastic 
scattering which can be used to study the optical 
properties of the detector and nuclear recoil. 

4. CCQE Events 

Charged current quasi-elastic interactions are 
fairly well measured in the MiniBooNE energy 
range on light targets. However, the cross section 
uncertainty on this process is an important error 
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contribution to oscillation searches, particularly 
for heavy target experiments like MiniBooNE. A 
summary of the existing data is shown in figure 
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Figure 2. Low energy CCQE cross section mea- 
surements vs. E v (GeV) [T3| . 

Quasi-elastic kinematics enable a precise deter- 
mination of the neutrino energy in v^(v e ) n — ► 
(i~(e~) p interactions. Neglecting corrections for 
the motion of the target nucleon, the neutrino en- 
ergy can be calculated from the measured energy 
and angle of the final state muon: 
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where M p is the proton mass, is the muon 
mass, Efj, is the muon energy, and 9^ is the muon 
angle with respect to the beam direction. The en- 
ergy resolution achievable by MiniBooNE assum- 
ing no non-CCQE background is ~10% at E v = 
1 GeV. 

5. Comparison of Monte Carlos 

Charged current quasi-elastic scattering is a 
simple process and therefore enables a straight 
forward comparison of cross section Monte Car- 
los. The NUANCE version 2 and version 3, 
NEUT, and NEUGEN Monte Carlos have been 



compared using the MiniBooNE flux, detector 
Monte Carlo, and reconstruction. This is the first 
such "apples-to-apples" comparison for CCQE in- 
teractions among cross section Monte Carlos. 

5.1. Theoretical Inputs 

The NUANCE, NEUT, and NEUGEN Monte 
Carlos have common theoretical inputs such as 
the LLewellyn-Smith free nucleon quasi-elastic 
cross section the Rein-Sehgal resonance cross 
section model 16 , and the standard deep inelas- 
tic scattering formula for high Q 2 ^7]. However 
there are non-trivial differences as well. These in- 
clude the implementation of the Fermi gas model 
for quasi-elastic interactions, the method for join- 
ing the resonance and deep inelastic scattering 
regions, and the treatment of final state interac- 
tions. 

The following is a brief summary of the salient 
theoretical inputs for these Monte Carlos. NU- 
ANCE version 2 uses dipole form factors, the 
Smith-Moniz Fermi gas model |18| . and mA 
1 .0 GeV/c 2 . NUANCE version 3 uses non-dipole 
form factors |19| . a new ir absorption model tuned 
on 7T data H3], and m A = 1.03 GeV/c 2 . NEUT 
uses dipole form factors, the Smith-Moniz Fermi 
gas model, and m A = 1.1 GeV/c 2 . NEUGEN uses 
dipole form factors, a tt absorption model tuned 
on v data, the Bodek-Ritchie modified Fermi gas 
model |21j. no nucleon re-scattering, and m A = 
1.032 GeV/c 2 . 

5.2. Normalization 

The predicted CCQE interaction rate is very 
similar between the NUANCE, NEUT, and NEU- 
GEN Monte Carlos. The CCQE percentage of 
the total number of events is 38.1% (NUANCE 
version 2), 39.8% (NUANCE version 3), 38.0% 
(NEUT), and 38.0% (NEUGEN). 

5.3. Kinematics 

The kinematic distributions for CCQE events 
considered here are reconstructed visible energy, 
reconstructed angle with respect to the neutrino 
beam direction, reconstructed quasi-elastic neu- 
trino energy, and reconstructed Q 2 . 

The visible energy is calculated from prompt 
light in the reconstruction, and is approximately 
equivalent to the muon kinetic energy. A com- 
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Figure 3. Left: reconstructed visible energy (GeV) of Monte Carlo CCQE events, normalized to unit 
area. Right: ratio of Monte Carlo predictions to NUANCE version 2 vs. reconstructed visible energy 
(GeV), shown with NUANCE version 2 CCQE cross section error band about 1.0. 



parison of the visible energy distributions from 
the various cross section Monte Carlos is shown 
in figure|3|(left), all normalized to unit area. NU- 
ANCE version 2 has a slightly harder spectrum, 
as can be seen in the ratio of the various Monte 
Carlos to the NUANCE version 2 prediction (fig- 
ure |3 right). The band about 1.0 indicates the 
CCQE cross section error currently assumed by 
MiniBooNE. 

The reconstructed angle of the final state muon 
with respect to the neutrino beam direction is 
shown in figure There are significant differ- 
ences between Monte Carlos in this distribution, 
particularly at high cos(6>) where NUANCE ver- 
sion 3, NEUT, and NEUGEN are all -5 - 10% 
lower than NUANCE version 2. 

The neutrino energy, calculated from the re- 
constructed muon energy and angle, is shown in 
figure For E$ E , both NEUT and NEUGEN 
exhibit softer spectra than NUANCE versions 2 
and 3. 

The reconstructed Q 2 distributions are shown 
in figure© The differences of - ±10% at Q 2 = 
between Monte Carlos are consistent with the 
comparisons of the cos((9) and visible energy dis- 
tributions. 

In general, it is striking that the NUANCE, 
NEUT, and NEUGEN Monte Carlos all agree 
within errors given the differences in their theo- 
retical inputs. The largest disagreements occur at 
high cos(#) and low Q 2 . It may be significant that 



this is the regime where the theoretically "messy" 
components of CCQE interactions, Pauli blocking 
and final state interactions, are most important. 

6. Comparison with MiniBooNE Data 

MiniBooNE data will eventually be used to 
measure the CCQE cross section, and even 
at this preliminary stage can be used to confront 
Monte Carlo predictions for event kinematics. To 
compare the cross section Monte Carlo predic- 
tions with MiniBooNE data, the same event se- 
lection cuts are applied to isolate CCQE events. 
The kinematic distributions discussed in the pre- 
vious section are then compared with the Mini- 
BooNE CCQE data set. 

6.1. CCQE Event Selection 

The CCQE selection criteria require that the 
event pass the cosmic background and fiducial 
volume cuts described in section [3 and that the 
event topology be consistent with a single muon- 
induced Cerenkov ring. The efficiency of the 
event selection is ~30%. For events generated 
within the 5 m fiducial volume the efficiency is 
~55%. The background is almost entirely due 
to charged current single pion production events, 

p — > fi~ p 7r + where the n + is produced below 
its Cerenkov threshold. The current MiniBooNE 
CCQE data set comprises ^60,000 events. 

The event selection performance varies slightly 
for the different Monte Carlos. A comparison of 
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Figure 4. Left: reconstructed cos(#) of Monte Carlo CCQE events, normalized to unit area. Right: 
ratio of Monte Carlo predictions to NUANCE version 2 vs. reconstructed cos(#), shown with NUANCE 
version 2 CCQE cross section error band about 1.0. 
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Figure 5. Left: reconstructed E~} E {GeV) of Monte Carlo CCQE events, normalized to unit area. 
Right: ratio of Monte Carlo predictions to NUANCE version 2 vs. reconstructed E® E (GeV), shown 
with NUANCE version 2 CCQE cross section error band about 1.0. 
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Figure 6. Left: reconstructed Q 2 (GeV 2 ) of Monte Carlo v„ CCQE events, normalized to unit area. 
Right: ratio of Monte Carlo predictions to NUANCE version 2 vs. reconstructed Q 2 , (GeV 2 ) shown with 
NUANCE version 2 CCQE cross section error band about 1.0. 
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the CCQE interaction rate before and after the 
event selection cuts is shown in table 1. The re- 
sulting CCQE purity and background contamina- 
tion is shown in table 2. 

Table 1 

Normalization comparison of Monte Carlos before 
and after CCQE event selection cuts. The ratio 
shown is the number of events with respect to 
NUANCE version 2. 

Monte ratio be- Efficiency ratio af- 

Carlo fore cuts of cuts % ter cuts 

NUANCE v2 LO 248 LO 

NUANCE v3 1.03 24.8 1.05 

NEUT 0.98 24.5 1.07 

NEUGEN 0.98 25.2 1.0 
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Figure 7. Reconstructed visible energy {GeV) of 
CCQE Monte Carlo and data, unit area nor- 
malized. 



Table 2 

CCQE purity comparison of Monte Carlos after 
event selection cuts. The dominant background 
is resonant single ir production, p — > fj," p it + . 

Monte CCQE Resonant % other 

Carlo % In % non-CCQE 

NUANCE v2 83 14 3 

NUANCE v3 80 16 4 

NEUT 78 13 9 

NEUGEN 80 16 4 



6.2. Kinematics Comparison 

The visible energy, cos(0), E® E , and Q 2 distri- 
butions after CCQE selection cuts are shown in 
figures [7| IS] |U and ^] respectively. The bands on 
the NUANCE version 2 Monte Carlo prediction 
indicate uncertainties from flux, CCQE cross sec- 
tion, and detector optical model variations. The 
data are shown with statistical errors. The NU- 
ANCE version 3, NEUT, and NEUGEN Monte 
Carlo predictions are overlaid. All distributions 
are normalized to unit area. 

The shape of the data in the visible energy 
distribution is significantly harder than the pre- 
dictions of the various cross section Monte Car- 
los, although the data are within the current er- 
ror band. The agreement between Monte Car- 
los and data is good for all cos(#) bins other 
than cos(0) ~ 1 where the NUANCE version 3, 
NEUT, and NEUGEN Monte Carlos agree with 



the data much better than NUANCE version 2. 
The data exhibit a somewhat harder spectrum 
in E® E than the Monte Carlo predictions. The 
most striking difference between data and Monte 
Carlo occurs near Q 2 ~ where the data are far 
lower than the predictions. 

It is interesting that although the cos(0) dis- 
tributions agree reasonably well for the data and 
some Monte Carlos, no Monte Carlo reproduces 
the low Q 2 behavior of the data. This may sug- 
gest a model deficiency common to all Monte Car- 
los considered here. 



7. Conclusions 

For the first time, "apples-to-apples" com- 
parisons have been made between cross section 
Monte Carlos. The NUANCE, NEUT, and NEU- 
GEN simulations agree surprisingly well with 
each other for CCQE interaction rates and kine- 
matics, given the differences in their theoretical 
inputs. There is disagreement at the 10% level at 
very low Q 2 between the MiniBooNE CCQE data 
set and all Monte Carlos considered here. In this 
regime, nuclear effects become most important. 
The differences between Monte Carlos and the 
MiniBooNE data may indicate a need for more 
sophisticated models of low Q 2 neutrino events. 
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Figure 8. Reconstructed cos(0) of CCQE 
Monte Carlo and data, unit area normalized. 
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Figure 9. Reconstructed E$ E (GeV) of ^ CCQE 
Monte Carlo and data, unit area normalized. 
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Figure 10. Reconstructed Q 2 (GeV 2 ) of 
CCQE Monte Carlo and data, unit area normal- 
ized. 
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